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Abstract 

Intersection angles of stable and unstable manifolds for area pre- 
serving mappings are numerically calculated by extremely accurate 



computation. With the use of multiprecision library the values of an- 



gle as small as 10 -400 are obtained. The singular dependence of the 
angle on the magnitude of hyperbolicity is confirmed. The power-law 
type prefactor with Stokes constant is also in good agreement with 
analytical estimation. 



1 Introduction 

Invariant manifolds are essential components in the phase space structure of 
Hamiltonian systems. A stable manifold W s and an unstable manifold W u 
of a fixed point (or a periodic point) is a set of points which asymptotically 
go to the fixed point as t — » oo and t — > — oo, respectively. These manifolds 
are invariant with time and are called invariant manifolds. 

Suppose we have an integrable Hamiltonian system Hq(I) where an unsta- 
ble manifold is smoothly connected to a stable manifold to form a separatrix. 
If we perturb the system as H = H (I) + eHi(I,(p) then the system is no 
longer integrable in general and the unstable manifold and stable manifold 
do not connect smoothly anymore. Rather they intersect transversally and 
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cross each other infinitely many times to produce a complex structure called 
homoclinic tangle, which eventually gives rise to chaos ||LL83| , |Tab89|| . Hence 
the intersection of invariant manifolds are quite important as the origin of 
Hamiltonian chaos. 




Figure 1: Intersection of unstable and stable manifolds and intersection angle 
9. The origin O is a hyperbolic fixed point. 



At the point of intersection between W u and W s we can define the "in- 
tersection angle" (which is also called as a splitting angle) 6 between these 
two manifolds. For several important cases it is estimated to be quite small 
and behave singularly with respect to the maginitude of hyperbolicity e of 
the fixed point |[LST89|| . Its leading order behavior for < e <C 1 is 



e 



C 2 

exp ( 



(i; 



where ci, c 2 (> 0) and //(> 0) are constants. For small values of e the ex- 
ponential factor makes the angle 9 extremely small. Recent asymptotic cal- 
culation of functional form of invariant manifolds also contain this singular 
factor [|HM93| , |TTJ94| , [TTJ98| , |NH96| , |HNK99f thus supporting the estima- 
tion. 

For area preserving mappings, the area enclosed by an unstable and a 
stable manifold is equal to the flux from inside to outside (or from outside 
to inside) of the separatrix per one time step. If the intersection angle gets 
small so does the area. Hence the smallness of the intersection angle is also 
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important in understanding the global dynamics and relaxation process in 
Hamiltonian systems. 

The singularly smallness of the intersection angle is, however, not yet 
confirmed thoroughly. Analytical estimation mentioned above is asymptotic 
approximation, and since it is asymptotic expression for e — > 0, we do not 
know from what value of e we can observe the singular behavior. 

Usually numerical computation can support analytical estimation. How- 
ever, ordinary numerical computation is useless in this case to check the 
singular behavior because of limited precision. Suppose we would like to 
confirm the proportionality of intersection angle as exp(— vr 2 /e). If e = 0.1 
the factor is about exp(— 100) ~ 10~ 43 . Let us write the normalized tangent 
vector and their components for stable and unstable manifolds as t s = (a, (3) 
and t u = (a + da, (3 + d(3) respectively, and the intersection angle as 9. Then 
| da | = 1 and \d(3\ = \a9\ up to the first order of 9. Hence if the angle is of 
the order of 10 -40 then the differences between values of the components of 
two vectors t s and t u are far beyond the range of ordinary double precision 
of 64-bit, where mantissa (fractional part) is at most 16 decimal digits. 

There are two methods which we may use for problems which require 
extremely high precision. One is "validated computation" ||N Y98|| and the 



other is "multiprecision" . "Validated computation" is a method to obtain 
rigorous results with numerical computation. It is based on operation for 
intervals rather than for numbers. Usual numerical computation on real 
numbers says e.g., x = 0.123, but this just mean that the value of x will be 
near 0.123 and it does not mean that the value of x is exactly 0.123. On 
the other hand, with the use of validated computation we can rigorously say 
that the value of x satisfies , for example, 0.122 < x < 0.124. 

Although validated computation is a powerful tool, it seems not to be 
all-purpose, and it seems to take time to reformulate the problem one want 
to solve to fit validated computation. Hence we adopt the second tool, mul- 
tiprecision library. 

In usual numerical computation the precision of variables are uniquely 
defined in each implementation of language (C , FORTRAN, etc.). On the other 
hand, using multiprecision library users can define the precision of variables 



arbitrarily. For various multiprecision tools see the hf loat web page [|Arn99 



Here we adopt a multiprecision library called cln |[Hai99|| which can be used 
in C++. 

The purpose of this paper is (i) to numerically confirm the singular de- 
pendence of the intersection angle using multiprecision library, (ii) to find 
the range of e where the singular dependence appears. 

There have been several attempts to numerically compute the singular 
behavior of the intersection angle. In the paper | |LST89| ] Lazutkin et al. 
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cites a numerical result by C. Simo as private communication and writes that 
the value of the angle by their analytical estimation is "in good accordance 
with numerical data". Also in | |A(JKRK92| |, Amick et al. computed the 
splitting distance of manifolds for a 4th order map x t+ 3 = x t + eg(x t+2 ), 
g(x) = x(l — x) with ordinary double precision. In this paper, we give a 
refined and concrete result for these pioneering work, by performing direct 
computation with multiprecision, 

This paper is organised as follows. In section 2 we introduce symplectic 
mappings for which we calculate the intersection angles. In section 3 we 
describe the method of computation, and the result is presented in section 
4. In section 5 we give summary and discussions. 



2 Models 

The first model we use is a "double-well map" ; 

0dw : (q,p) i-> (q',p') , p' = p- e(2q 3 - q) , q' = q + ep . 

This map has a hyperbolic fixed point at (0,0). 

For this mapping the intersection angle is estimated to be [ |NH96 | 



n C dw / 7T 

#dw = — r- exp 



5.00 x 10 4 , 



after change of coordinate^. 

The second model is the well-known standard map ||Chi79 



K 

frst : (q, P) >-> (q, P) , P = P + 7T sin (2vrg) , q =q + p . 

Z7T 



(2) 



(3) 



(4) 



This map has hyperbolic fixed points at (n, 0) (n : integer). For standard 
map the intersection angle is estimated to be ||LST89|| 



9 st = — exp 
K 



7T 



c st = 1118.827706.. 



(5) 



3 Method of calculation 

Here we describe the method of calculation of intersection angle for the double 
well map (fj). The method for standard map is similar. 



1 The paper [ |NH96 ] adopts different notation for the map (0) hence the value of the 
angle 9 and the prefactor Cd w is rescaled appropreately in eq. (pf). 
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In short, we choose two approximate points P± for an intersection point, 
and approximate the intersection angle by the ones between the two mani- 
folds which are on P±. 



3.1 finding an interval which contains PIP 

First we search for an intersection point of W u and W s . Following an argu- 
ment similar to [|GLS94 ], it is easily shown that one of the intersection points 
called principal intersection point (PIP) is on the g-axis. 




Figure 2: True stable and unstable manifolds (dashed lines) and approximate 
manifolds (solid lines) and tangent vectors 



Suppose we choose an initial point on g-axis and update the point back- 
ward in time. If the point is near but not on the stable manifold it will 
approach to the hyperbolic fixed point and then goes away from it. The 
direction the point goes away shows whether the initial point is on the right 
or left of the true intersection point. Hence we can set an interval on q- 
axis which include the intersection point between the stable and unstable 
manifolds. 

Then we make the interval narrower by bisection. We call the final interval 
as (g_, q + ). In actual calculation we set the final width of the interval as 

\ q+ -q_\ = 1(T 100 for e > 0.02 (6) 
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and 



\ q+ - q _\ = 1(T 8U0 for e = 0.01 (7) 

We will use the points (g_, 0) and (q + , 0) as approximations for the intersec- 
tion point and denote as P_ and P + respectively, i.e., 

P± = (g±,0) (8) 

3.2 construction of tangent vectors in the linear region 

Now we are going to make tangent vectors at these approximate points P±. 

Although we do not know directly the tangent vectors at P±, we can 
calculate them from the linear neighbourhood around the hyperbolic fixed 
point (0,0). 

First we update the points P± forward or backward in time; 

p^ = ^(p±), Pi = <t>t( p ±) (9) 

where the time steps N u and N s are sufficiently large for |P"| <C 1 and 
\P±\ 1. In actual calculation we determined N u and N s so that P± and 
P± be the points which are nearest to the hyperbolic fixed point (0,0) on the 
orbits (P±) and <p^ w (P±), respectively. 

Since the points P± and P± are quite near to the fixed point (0,0) the 
map 0dw can be approximated by its linearized map; 

dw (g,p) =D0 dw (O,O) r) + 0(q 2 ,p 2 ) . (10) 

Hence the tangent vector of the manifolds which passes P u or P s can be 
approximated by the invariant curve of the linearized map D(j)(0, 0), that is, 
hyperbola. If we denote 



c — b a — d 
Then the curve 

$(q,p) = (b - c)q 2 + 2dqp + (b + c)p 2 = const. (12) 

is invariant under dw ||Grc68|| . Hence we can obtain the tangent vector 
and If of unstable and stable manifolds as the tangent vector of the 
hyperbola (|12"D 

= «P® -Qq- ) (Pl)n = uors. (13) 
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3.3 calculation of intersection angles 



The tangent vectors at approximate intersection point P± is obtained by 
pulling back the tangent vectors t u , t s to the points P± as 



5* ± 



(14) 



Approximation values for intersection angle 9 at each approximate inter- 
section point P± is obtained as 



COS^H 



v u ± ■ V ± 



(15) 



4 Results 

Now we present the results. Computation is performed on a Pentium-II pro- 
cessor PC with Linux operating system (Vine Linux 1.0b) with GNU C/C++ 
compiler g++ 2.7.2.3 and egcs 1.0.3-release, and the version of multiprecision 
library is cln 1.0.1. The computation time is typically several hundred times 
longer than ordinary numerical computation with double precision. 

The result for double well map @ is summarised in the Table [l] and 
Fig|| Calculation is performed with 1000 decimal digits for e > 0.02 and 
4000 digits for e = 0.01. 



e 


9_/n 


9 + /7l 


0.6 


0.00620123 


0.00620123 


0.4 


2.03305e-05 


2.03305e-05 


0.2 


1.61224e-14 


1.61224e-14 


0.1 


2.0686e-34 


2.0686e-34 


0.08 


1.23026e-44 


1.23026e-44 


0.06 


7.24697e-62 


7.24697e-62 


0.04 


1.06168e-96 


1.06168e-96 


0.02 


2.58525e-202 


2.38495e-202 


0.01 


1.16442e-414 


1.16442e-414 



Table 1: Intersection angle 9 for double well map. 1000 decimal digits for 
e > 0.02 and 4000 decimal digits for e = 0.01. \q L -q R \ = 10" 100 for s > 0.02 
and 10" 800 for e = 0.01 . 



Here we can see that our numerical results agree quite well with the 
estimation (eq.(H)). It is clear that the data cannot be fitted by power law. 
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intersection angle for double-well map 




Figure 3: Intersection angle for double well map (fj). Values are in Table(l|. 
The dotted line represents analytical estimation for the singular behavior 
exp(-7r 2 /£)(eq.(§)). 



Also we note that the singular have behavior appears for the value of e as 
large as e — 0.1. 

The difference of the values between 6L and 8 + is an indicator of error. 
Table [I] shows the error is small. 

Next we see the results for standard map (|]). For standard map we have 
the principal intersection point on the line q = it [|GLS94|| , hence we take its 



approximation points as P± = (n,p±). Calculation is performed with 1000 
decimal digits. For standard map in the notation (^) the parameter e is equal 

to Vk. 

From the figure [| again we see that the intersection angles agree quite 
well with analytical estimation (|5|). Note that we apparently have agreement 
even above the K c , where the last KAM torus is broken and the phase space 
has global chaotic sea. 

Finally we show the prefactors of the singular behavior {ciE~^ in eq. ([!]).) 
Here we examine whether the prefactor is in fact a power law type with 
previously estimated exponent [i. Also the value of the constant c\ is in- 
teresting to check, because c\ contains a constant called Stokes constant, 
which represents the change in the form of asymptotic expansion of the 
manifolds [TTJ91 , [TTJ98| , |HNK99| , [RH96 ] 



Fig. |^ shows the e dependence of the value 9 exp(7r 2 /e). Here we see a 
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K 


6-/7T 


6 + /tt 


1.5 


0.0692401 


0.0692401 


1.0 


0.0250374 


0.0250374 


0.8 


0.0114877 


0.0114877 


0.4 


0.000324148 


0.000324148 


0.2, 


l.lb4U/e-Ub 


l.lo4lHe-Uo 


0.1 


2.6934e-10 


2.6934e-10 


0.08 


8.60911e-12 


8.60911e-12 


0.04 


9.47743e-18 


9.47743e-18 


0.02 


2.58424e-26 


2.58424e-26 


0.01 


1.46928e-38 


1.46928e-38 


0.008 


1.60962e-43 


1.60962e-43 


0.004 


4.59792e-63 


4.59792e-63 


0.002 


7.84703e-91 


7.84703e-91 


0.001 


3.14686e-130 


3.14686e-130 



Table 2: Intersection angle for standard map. 1000 decimal digits are used. 
\ p+ - p _\ = io- 100 . 



— 1e-60 
Q. 



c 



Intersection angle for standard map 



"st2.data" using 1 :2 — I— 
exp(-pi**2/sqrt(x)) — 1— 



K 



Figure 4: Intersection angle for standard map. Values are in Table |2|. 
The dotted line represents analytical estimation for the singular behavior 
ex V (-n 2 /VK) (eq.(D). 
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prefactor of splitting angle for double-well map 



"data9.data" using 1 :4 — I— 
41 849. 5"eps"(-5. 04935) 

5.00*1 0"4-eps**(-5) 



eps 



Figure 5: Prefactor of intersection angle for double-well map (|2]). Dashed line 
represents a least square fitting of the data for 0.01 < e < 0.1, and dotted 
line represents an analytical estimation Cd-w/e 5 (|3|). 



prefactor of intersectino angle for standard map 

"st2.data" using 1 :4 — |— 
924.551/K"(1. 02875) — 

1118.8277/K 



-•4=.. 



C\J 100000 

<_ 
a. 

CL 

X 

10000 

jd 

cB 1000 [ 



+ 



Figure 6: Prefactor of intersection angle of standard map. Dashed line rep- 
resents a least square fitting of the data for 0.001 < K < 0.1 , and dotted 
line represents an analytical estimation c st /K (^). 
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clear power-law dependence for small values of e. A least square fitting for 
0.01 < e < 0.1 gives 

9 ■ exp(7r 2 /e) = 41894.5/ £ 5 - 04935 , (16) 

which is in good agreement with the estimation (^). The difference in the 
exponent ( 5.04935 for our numerical result and 5 for |[NH96|| ) is negligible. 
Also the prefactor 41894.5 is in good agreement with the estimated value of 
5.00 x 10 4 . 

Fig. |^ is for the standard map, where we see a fitting for 0.001 < K < 0.1 

6/7texp(7t 2 /VK) = 924.551/K 1 ' 02875 , (17) 

again in good agreement with the estimation The difference between our 
numerical result and the estimation is negligible, as shown in the figure. 



5 Summary and discussion 

In this paper we have numerically computed the intersection angles between 
stable and unstable manifolds for two area preserving mappings. Computa- 
tion is performed with extremely high precision ( with 1000 or 4000 decimal 
digits) with the use of multiprecision library cln to obtain the angle as small 
as 10~ 400 . The singular behavior previously obtained by analytical estimation 
is confirmed. Hence we have confirmed that the singular small behavior of 
the intersection angle 9 cx exp(— c/e) indeed exists for area preserving map- 
pings. Thus the flux in the phase space is quite small for small values of e, 
which will make the relaxation process quite slow. In addition, although the 
analytical estimation is obtained as the asymptotic approximation formula 
for e — > , our numerical results show that the singular dependence appears 
from rather large value of e. 

Also we have measured the prefactor of the singular factor, and confirmed 
that the analytical estimations for prefactors are correct. Exponent of the 
power law is in almost perfect agreement. It appears constants in the pref- 
actors slightly differ from the ones in analytical estimation. This does not 
necessarily imply that the estimation need to be refined. For the analytical 
estimation contains a free parameter called switching function 5(0) which 



take a value between and 1 ||NH96|| . Although the estimation is given for 
5(0) = 1/2, 5(0) can be some other value. In fact ||NH96|| takes the value 
5(0) = 0.3 when comparing their functional form of unstable manifold with 
numerically obtained orbit. Note that within the range of parameters com- 
puted the constant factor do not appear to depend on e or K. 
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The singularly smallness of intersection angle was one of several things in 
the study of Hamiltonian chaos which are predicted but not yet confirmed 
precisely. Among others are Nekhoroshev bound ||Nek77| , [LJN 92| | and Arnold 
diffusion ||Arn64 |. In both cases singular factor appears and it is believed that 
the singular behavior accounts for slow motion in high-dimensional nearly in- 
tegrable Hamiltonian systems ||KK90|| . Extreme precision computation using 
multiprecision library will reveal some new properties about the slow dynam- 
ics of high-dimensional Hamiltonian systems. 
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